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We review current ideas on entropy production during the different stages of a relativis- 
tic nuclear collision. This includes recent results on decoherence entropy and the entropy 
produced during the hydrodynamic phase by viscous effects. We start by a discussion 
of decoherence caused by gluon bremsstrahlung in the very first interactions of gluons 
from the colliding nuclei. We then present a general framework, based on the Husimi 
distribution function, for the calculation of entropy growth in quantum field theories, 
which is applicable to the early ("glasma") phase of the collision during which most of 
the entropy is generated. The entropy calculated from the Husimi distribution exhibits 
linear growth when the quantum field contains unstable modes and the growth rate is 
asymptotically equal to the Kolmogorov- Sinai' (KS) entropy. We outline how the ap- 
proach can be used to investigate the problem of entropy production in a relativistic 
heavy-ion reaction from first principles. We show that the same result can be obtained 
in the framework of a completely different approach called eigenstate thermalization 
hypothesis. Finally we discuss some recent results on entropy production in the strong 
coupling limit, as obtained from AdS/CFT duality. 

1. Overview 

1.1. The General Problem 

High-energy heavy-ion coUisions show a vary rapid transition from a quantum me- 
chanical initial state to a thermalized state known as a quark-gluon plasma. This 
poses a fundamental, conceptual theory problem, namely how the necessary produc- 
tion of entropy can be reconciled with the T-invariance of quantum chromo dynamics 
(QCD). One could argue in analogy to the situation in the Einstein- Rosen-Podolsky 
paradox that all produced particles stay in a highly entangled state which only 
collapses when a measurement is performed and that entropy production occurs 
only in that moment. Nevertheless, a hydrodynamic description of the quark-gluon 
plasma seems to be possible long before any measurement takes place. The reso- 
lution of this apparent contradictionis of interest for a wide spectrum of physical 
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phenomena, e. g., reheating of the early universe after the end of the presumed pe- 
riod of inflationary growth, when the universe transitioned from an exponentially 
inflated vacuum bubble to the hot cosmos we observe today. Another example is 
the gravitational collapse of matter into a black hole, which behaves as if it were 
a thermal object. In this review we will focus on the thermalization of strongly in- 
teracting matter after the collision of two relativistic heavy nuclei, thereby forming 
a quark-gluon plasma, a process which has the immense advantage that it can be 
studied in great detail in laboratory experiments. For example it is known that the 
quark gluon matter in a heavy-ion collision behaves like a nearly ideal thermal fluid 
already after 1 — 2 fm/c. 

Justifying the use of thermodynamics in the context of quantum field theory is 
an extremely difficult task in any case, even for phenomena described by quantum 
electrodynamics (QED), where the quantum field theory is under excellent control. 
For QCD the situation is even more difficult due to the strong coupling dynamics, 
which permits analytical calculations only in very much constrained situations. To 
illustrate the extent of the problem let us note that even in quantum mechanics 
the problem of thermalization is far from being settled. A nice example is a recent 
study^, where a finite one-dimensional quantum system was shown to thermalize 
only when the so called "eigenstate thermalization hypothesis" (ETH^) is satis- 
fied. Rigol^ stresses that quantum and classical thermal states have very different 
natures. ETH claims that each eigenstate of the Hamiltonian is related to a ther- 
mal state which is reached by the "de-phasing" of the coherent components of the 
quantum state. 

While we do not want to discuss the merits of this specific concept in more 
detail, it serves as an example illustrating the fact that the de-phasing or decoher- 
ence of a quantum state constitutes one way (if not the way) to generate entropy. 
Generated entropy generally is proportional to the number of lost bits of informa- 
tion. Consequently phases carry an important part of the information contained in 
a quantum state, their loss generates a significant amount of entropy. 

A general mechanism for losing information about the state of a quantum sys- 
tem is the entanglement of its wave function with that of its environment. This 
mechanism obviously does not apply to an isolated system. The conservation of the 
von Neumann entropy for an isolated system implies that any information encoded 
in its quantum state is never lost, at least in principle. However, the information 
may become encoded in observables of such complexity or requiring such a degree 
of experimental precision that its retrieval is impractical. A standard approach to 
describe this effective loss of information in an isolated system is coarse graining. 
By averaging over finite regions of phase space, information is lost, and thus entropy 
is produced. 

Often, coarse graining is discussed in connection with measurements performed 
on a system. However, in Section 2 we will discuss how even in the absence of 
physical measurements the uncertainty principle implies a natural, measurement 
independent coarse graining. It thus implies effective entropy production for non- 
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linear, time-reversal invariant systems, in which individual trajectories approach 
all points in some region of phase space, while all ensembles of trajectories respect 
Liouville's theorem, i.e. occupy a constant phase space volume. 

This argument also applies to the information residing in the phases between the 
different components of a quantum system, which are smeared by coarse graining. 
Thus decoherence or de-phasing plays a crucial role for the discussion of ther- 
malization and entropy production, which naturally leads to the question whether 
decoherence alone is sufficient to reach thermalization. This is not the case, because 
decoherence may lead to systems with statistical features, which have energy and 
momentum distributions that are far from equilibrium. 



1.2. Relativistic Heavy Ion Collisions 

Heavy-ion collisions provide a good example for such systems. After decoherence 
of the initial parton distributions in the colliding nuclei one is left with some sort 
of quark-gluon matter that can, for most part, be treated statistically and semi- 
classically. The natural question is whether the equations of fluid dynamics can be 
applied to describe the subsequent evolution of this matter. Fluid dynamics is an 
effective theory for the time evolution of the expectation value of the components 
of the energy- momentum tensor T^j^, which encodes the dynamics of the system 
in terms of macroscopic quantities: energy density, pressure, collective flow veloc- 
ity, and the equation of state. Its applicability requires that spatial and temporal 
gradients of these quantities are small, and that the deviations of the microscopic 
statistical state of the system from local thermal equilibrium are small and can be 
represented by viscous (gradient) corrections describing anisotropics of the stress 
tensor. 

The early-time energy momentum tensor for heavy-ion collisions (colliding along 
the z axis) obtained by solving the classical Yang-Mills equations of motion has the 
approximate form^ 
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with negative longitudinal pressure, which makes it unsuitable for hydro dynamical 
calculations. The negative longitudinal pressure is an expression of the fact that 
Gauss' law does not permit all components of the gauge field to decohere, because 
they are not independent of the field generating color charges. In this case, sec- 
ondary gluon production depletes the remaining coherent longitudinal fields and 
eventually generates positive longitudinal pressure. Only then can the system es- 
tablish acceptable initial conditions for the following hydrodynamical evolution. 
The particle production process and the interaction among the produced particles 
therefore must be included in the description of the thermalization process. 
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A natural question to ask in this situation is, how quickly the components of the 
stress tensor develop a sufficient degree of isotropy in the locally co-moving frame 
for a hydrodynamical description to be valid, and how this isotropization time^ riso, 
compares to the decoherence time, Tdec- Note that isotropization is a less strin- 
gent requirement than thermal equilibration. In principle, an isotropic stress tensor 
could be generated by an isotropic, yet strongly non-thermal gluon distribution. 
Indeed, if one considers a dilute ensemble of de-phased particles, which is described 
by the Boltzmann equation, in a state near thermal equilibrium, its relaxation to 
equilibrium is governed by an infinite set of relaxation times and their associated 
modes. Because of the local conservation of energy and momentum, the viscous 
modes are usually those which survive at late times, and thus entropy growth at 
late time is usually controlled by the viscosities and the thermal conductivity. But 
this does not have to be the case when the system is far off equilibrium, where 
unstable, exponentially growing modes may exist that do not permit a linearized 
description. Such processes can somtimes be interpreted as run-away (stimulated) 
particle production. 

While particle production is usually not interpreted as decoherence, it may in 
fact be related to it. For example, bremsstrahlung can be understood as decoher- 
ence of the coherent, quasi-real Weizsacker- Williams fields associated with any fast 
moving charge. While such radiation does not result in exponential growth in QED 
in the vacuum, it can do so in the presence of strong fields (e.g. in a free electron 
laser) or in the nonlinear dynamics of an electromagnetic plasma or a quark-gluon 
plasma far off equilibrium. The interpretation of how entropy is generated in such 
processes may then depend on the theoretical and conceptual framework in which 
the processes that lead to thermalization are described. 

1.3. General Considerations 

The question of whether, how, and how fast a highly excited state of QCD matter 
equilibrates thermally is thus a multi- faceted one, and we certainly do not profess 
to fully understand all relevant points of view. Below, we will outline and discuss 
those we are aware of, without claiming completeness. A general framework for 
the study of thermalization processes starts from the observation that thermal 
equilibrium corresponds to the state of highest entropy. Our approach thus focuses 
on the rate at which entropy is produced in the system. As already mentioned, the 
von Neumann entropy of an isolated quantum system never increases, because the 
time evolution of an isolated quantum state is unitary. However, when we speak of 
the thermalization of an isolated system, we mean that its state "looks" thermal to 
most practically feasible experiments that measure a select number of observables. 
Owing to the lack of complete observation we replace the microscopic, possibly pure 
quantum state of the system by a density matrix that traces over all unmeasured 
observables. By this coarse graining procedure we assign an entropy to the system, 
which can grow as a function of time if the microscopic structure of the quantum 
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state gets hidden in more and more complex, and hence practicahy unobservable, 
degrees of freedom. The coarse grained, apparent entropy is then a measure of the 
complexity of the microscopic state of the system. 




Fig. 1. Sketch of entropy growth in a generic non-hnear system. Initiahy the behavior is non- 
universal, depending on the specific initial conditions, then one encounters a phase of linear entropy 
growth and finally an asymptotic approach of complete thermal equilibrium 



The generic time evolution of the coarse grained entropy for a dynamical system 
with ergodic properties is sketched in Fig{l] At very early times, the behavior is 
extremely sensitive to details of the initial state and will depend on the initial 
phase correlations and the occupation probabilities of the various eigenmodes of 
the system. If all modes are randomly occupied the length of this period can be 
estimated^ to be of order a/TV, where N is the number of independent modes. The 
period of approximately linear rise of the coarse grained entropy is governed by 
the sum of all positive Lyapunov exponents of the system in the classical limit, 
the so-called Kolmogorov- Sinai' (KS) entropy growth rate. The rise starts when the 
amplitude of the unstable modes has outgrown the average mode occupation, and 
it terminates when the coarse grained entropy approaches the equilibrium entropy. 
While the rate of growth of the entropy is an intrinsic property of the system, the 
total duration of the quasi-linear period depends on the specific nature of the initial 
state. Finally, the last phase is governed by the relaxation times of small deviations 
from equilibrium. 

The dynamical evolution of the quark-gluon plasma created in a relativistic 
heavy ion collision is commonly described by relativistic viscous hydrodynamics, 
which is well suited to describe the last, asymptotic phase sketched in Fig{l] Ob- 
viously, the question, when the transition from the microscopic description to the 
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hydrodynamical one should be taken, has no unique answer but has to be decided 
on the basis of heuristic and practical considerations. In principle, the initial con- 
ditions for the hydrodynamic simulation will depend on which moment is chosen 
as the starting time. However, if hydrodynamics applies at all, the late time results 
of the simulation will depend only little on when exactly the transition from a mi- 
croscopic description to hydrodynamics is made. According to our discussion, the 
only constraint is that this time should be chosen close to the end of the phase of 
linear entropy growth. 

In fact, the situation is even more involved. In practice, the information about 
the supposedly thermal character of the final state is obtained from single-particle 
observables, including emitted particle yields, spectra, and angular distributions, 
as well as of selected two-particle correlation functions, such as balance functions. 
These observables are all of hadronic nature; thus, one also has to take into ac- 
count the entropy production which occurs during hadronization, as well as that 
caused by interactions in the expanding hadron gas before freeze-out. These final 
one-particle or limited two-particle observables trace over a very large number of 
degrees of freedom and thus imply a highly coarse grained description in terms of 
the fundamental degrees of freedom. However, we are not so much interested in this 
late state as in the question when the transition to a hydrodynamical description 
is justified early in the reaction. For this question, a coarse grained entropy defini- 
tion that is dictated by the exigencies of final-state measurements, long after the 
collision is over, makes little sense, because a thus defined entropy may have very 
little relevance for the dynamics governing the time evolution of the system. This 
suggests that it may be useful to study entropy growth under a minimal amount 
of coarse graining (as imposed by the uncertainty principle, see below) and to use 
this approach to estimate the rate at which the coarse grained entropy approaches 
its thermal equilibrium value. 

The different conceptual approaches to the thermalization problem are to 
some extent complementary. Obviously, when the system truly equilibrates, all ap- 
proaches will indicate thermalization. However, isotropy of its macroscopic proper- 
ties alone is a poor measure of thermal equilibrium of a quantum system, nor does 
the growth of the coarse grained entropy ensure that the system approaches local 
isotropy. A fiuid dynamical description with a thermal equation of state requires 
both. In this review, we focus on the progress of our understanding of how and 
when the apparent entropy in the final state of a relativistic heavy ion collision 
is generated by the internal dynamics of the highly excited QCD matter. How- 
ever, we caution the reader that, while substantial progress has been made over 
the last years, we are still far from having a complete answer to the problem of 
thermalization. 
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2. Entropic History of a Relativistic Heavy Ion Collisions 
2.1. Stages of entropy production 

As already mentioned, much of our information on the behavior and properties of 
hot QCD matter is derived from measurements of the particle yields and spectra in 
the final state of relativistic heavy ion collisions and their interpretation in terms 
of thermodynamic and hydrodynamic concepts. At the energies available at the 
Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) one 
of the most relevant quantities is the azimuthal quadrupole anisotropy of the collec- 
tive flow, usually called elliptic flow and denoted as V2. An excellent agreement of 
hydrodynamical calculations of the flow anisotropy of the matter produced in the 
nuclear collisions with the elliptic flow measurements requires the assumption of a 
rapid thermal equilibration of the matter on a time scale of the order of 1 fm/c. 
An important problem in the description of relativistic heavy ion reactions is thus 
to understand how the produced matter equilibrates so quickly, i.e., when and how 
entropy is created in the reaction. In principle, one can distinguish five different 
stages of entropy production: 

(1) Decoherence of the initial nuclear wave functions; 

(2) Thermalization of the partonic plasma ("glasma"); 

(3) Dissipation due to shear viscosity in the hydrodynamical expansion; 

(4) Hadronization accompanied by large bulk viscosity; 

(5) Viscous hadronic freeze-out. 



Decoherence 
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Fig. 2. Entropic history of a central Au+Au collision at top RHIC energy. The values of dS/dy 
indicate the entropy per unit rapidity reached at the end of various collision stages based on 
experimental data and model estimates. 
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Our present incomplete knowledge about the contribution of these different 
stages to the final entropy is depicted in Fig. [2] Before we discuss these stages in 
turn, we note that the final entropy per unit rapidity dS/dy is one of the best 
known quantities in relativist ic heavy ion physics. In the case of Au+Au collisions 
at RHIC the entropy distribution dS/dy at particle freeze-out can be determined 
from an analysis of the final hadron spectra in combination with the information on 
the source radius derived from identical particle (HBT) correlations.'^ The slightly 
extrapolated result for the 6% most central Au+Au collisions at y^^NN = 200 
GeV is {dS/dy)i = 5600 ± 500 at midrapidity.l^ Alternatively, the final entropy can 
be deduced from the measured hadron abundances, combined with the calculated 
entropy per particle for a hadron gas in chemical equilibrium at Tc ^ 160 MeV, 
S/N ^ 7.25, which yields the result^ {dS/dy)ch = 5100 ± 200. The 10 percent 
difference between these values can be attributed to the entropy production during 
the hadronic freeze-out and reflects the entropy production due to the decay of 
excited hadronic resonance states as well as the contribution from the substantial 
shear viscosity of a thermal hadron gas.^^ 



2.2. Decoherence 

Entropy production by decoherence dominates the very first stage of a heavy ion 
collision. The loss of coherence is measured by the decay of off-diagonal elements 
of the density matrix p describing the system. A practical way of investigating this 
process is to calculate the decay rate of the quantity Trp^/[Trp]^. The time scale 
of the decoherence of the initial nuclear wave function has been studied by this 
method in the color glass condensate model (CGC).^ In this model, the initial 
nuclear gluon distributions are coherent over a wide range of longitudinal momenta 
kL (in the beam direction) due to the high degree of Lorentz contraction of the 
fast moving nuclei. The elementary process driving the decay of the off-diagonal 
elements of the density matrix is the fusion of one gluon from each nucleus, which 
removes gluons from the initial-state wave functions. The relevant time scale in this 
process is set by the saturation scale Qs of the nuclear gluon distribution. Qs is a 
measure of the average transverse momentum of gluons in a fast moving nucleus. 
Gluons with transverse momentum kr < Qs ^ ^ — GeV can be thought of as 
components of a quasi-classical gauge field. A careful analysis yields the result that 
the characteristic decoherence time is r^ec = cQ^^, where c denotes a calculable 
constant close to unity.^ While this result is expected on dimensional grounds, 
the fact that c = 0(1) is important because it shows that decoherence occurs over 
a time of less than 0.2 fm/c. 

Here we do not want to present the complete chain of arguments leading to 
this result, but only provide a few details showing how the CGC model enters into 
the calculation. We study perturbative gluon interactions between the gluons in a 
nucleus at rest (nucleus 1) and a fast moving one (nucleus 2). The time evolution 
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of the density matrix describing nucleus 1 reads'^ 

PAjt) = Y.f didt' ifjy (t>^,B(0)^f^,T(^') (2) 

B,B ^ 

The CGC model enters by specifying the correlator {A^{p)Ap^{q))2 between gauge 
fields in nucleus 2 as 

/ 2 P Q P_\_ 

where 

F{p±) = [ d^xe-'^---^^^^^ (1 - e-9'^c/(8^)M^x^ini/(|x|A)\ _ 

This form implies decoherence in transverse momentum but coherence in longitu- 
dinal momentum. For the gluon field correlator in nucleus 1 (at rest) one makes a 
simple Gaussian ansatz 

=5-,Ma{2^f5\k-mee-^'^''''+'''^ . (5) 



It turns out that a finite result is only obtained if the following running of the QCD 
coupling constant is usec^ 

/ gHA')gHl/x') , (6) 

where A is the infrared cutoff of the CGC, which cancels in the ratio Trp^/[Trp]^. 
This ratio then turns out to be inversely proportional to the life-time of the system, 
and one can define Tdec as the time after which this ratio has dropped by a factor 
l/e. 

The amount of entropy created from the loss of coherence of the components of 
the nuclear gluon wave functions can be estimated as follows. The complete deco- 
herence of a coherent quantum state with average occupancy n results in a mixed 
state with entropy^^ 5'dec ^ (ln27rn + l)/2. The decoherence occurs in transverse 
domains of spatial size Q~^, which are causally disconnected during the decoher- 
ence process. The number of causally disconnected transverse domains in a central 
Au+Au collision is of order (Qs^)"^, where R denotes the nuclear radius. Ac- 
counting for a longitudinal coherence length Ay ^ 1/as and using n ~ l/o^s one 
findsIMH 

(dS/dy)^,, ^ ^QlR^as{\n27r/as + 1) ^ 1,500 (7) 



at midrapidity for a central Au+Au collision. This is roughly one quarter of the 
observed final entropy per unit rapidity. 
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2.3. Viscous flow 

The other stage which is reasonably weh understood is the hydrodynamic expan- 
sion stage. RHIC data on the ehiptic flow of hadrons in non-central collisions 
indicate that the shear viscosity during this phase is small. The bounds on the 
shear viscosity rj compatible with the RHIC data lie characteristically in the range 
1 < Atttj/s < 2.5F^I"'^^I"'^^I"'^^I"'^^I"'^^] Such a small viscosity cannot increase the total 
entropy of the expanding fluid by much. Systematic studies of the entropy growth 
in longitudinally boost-invariant, viscous fluid dynamics for different initial con- 
ditions were performed by Fries et al^ The viscous stress-energy tensor in the 
locally co-moving frame has the form T^^ = diag(£:. Pi, Pi, P — z). The transverse 
and longitudinal pressures can be decomposed into the thermodynamic pressure P 
and the components ^ and 11 of the shear and bulk stress: P± = P + 11 + 

p^ = p + n - 

At early times, ^ and II are given by the initial conditions for the stress ten- 
sor established during the thermalization process. In view of the dilution effect of 
the longitudinal expansion imprinted on the matter by the collision, the pressure 
components are expected to satisfy the ordering P^ < P±. The smallest physically 
meaningful value of the longitudinal pressure is P^ = —Pi, which arises when the 
matter is completely in the form of coherent longitudinal fields 1^1221 However, the 
fluid dynamical approach does not apply to such a situation, because it assumes 
small deviations from local equilibrium. The minimal condition for fluid dynamics to 
be applicable is the mechanical stability condition: P^, Pi > 0, but the regime where 
the fluid dynamical description can be trusted only begins when ^/P, H/P<:1. 

The equations governing the longitudinal expansion of the medium areEHMHEH 

|£ = -i(£ + P + n-$), (8) 

or T 

^dT 3r ^ ' St ^ 2r]^ ' 

rn^ = -^-n(r). (10) 

rj and ( are the shear and bulk viscosity, respectively. The relaxation times for 
the shear and bulk stress, r^^ and rn, were assumed to be those derived from the 
Boltzmann equation, A was taken from the supersymmetric gauge theory (see Fries 
et al^foT details). The entropy per unit rapidity and transverse area, dS/ (dydA) = 
rs obeys the equation:^ 



d{rs) _ T /3^2 , H 



2 ^ 

(11) 



dr T \4r] C 

Fries et al^ started the hydrodynamical simulations after an initial time tq = 
0.3 fm/c, which is compatible with the expected decoherence time but considerably 
shorter than the equilibration times estimated from the elliptic flow analysis. The 
initial energy density was fixed at Sq = s{ro) = 50 GeV/fm^. The early starting 
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Fig. 3. (a) Relative transverse and longitudinal pressure, P^/P (dashed) and Pz/P (solid), as 
functions of time r for the initial conditions (i) (black), (ii) (blue), and (iii) (red), (b) Relative 
entropy production from bulk and shear stress, Sn/Sf (solid) and S^/Sf (dashed) for the same 
scenarios. For further details, see Fries et ah— 



time was chosen to explore, in a schematic way, how long it would take to reach the 
domain of applicability of the hydrodynamical description defined by the conditions 

- ^' 1 \Pz-P\ 1 

< o; < o- (12) 



P - 2' P - 2 

Fries et al^^ considered three different initial conditions: 

(i) n(ro) = ^(to) = 0, corresponding to thermal equilibrium at tq; 

(ii) n(ro) = —C,{Tq)/tq and ^{tq) = 47^(To)/(3ro), corresponding to the shear stress 
of the relativistic Navier-Stokes theory; 

(iii) n(ro) = —C{To)/ro as in (ii), but ^{tq) = P{tq) + n(ro) for vanishing initial 
longitudinal pressure. 

Typical results from these calculations are shown in Fig [s] The lower panel (a) 



indicates that the condition (12) is approximately satisfied for r > Teq ~ 1 fm/c 
independently of the choice of the initial condition for the stress tensor. On the other 
hand, the upper panel (b) shows that most of the viscous entropy production occurs 
before this time during a period when the validity of the hydrodynamical approach is 
questionable. The relative contribution to the final entropy from viscous forces after 
Teq does uot cxcccd 10 percent for the parameter choices explored in this calculation. 
The same conclusion was reached by Song and Heinz^ in a two-dimensional boost- 
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Fig. 4. Relative entropy increase due to bulk and shear stress, AS/Sq, for two-dimensional boost 
invariant viscous hydrodynamical expansion of a quark-gluon plasma for three different equations 
of state. The entropy increases by approximately 10 percent, if the shear visscosity is equal to the 
KSS bound r]/s = I/Att. For further details, see Song and Heinz^-*^. 



invariant calculation (see Fig. [4|. We may thus conclude conservatively that the 
entropy at the moment of local equilibration must be {dS/dy)^^^ > 4,500. These 
considerations tell us that approximately half of the final entropy must be generated 
during the thermalization process, which cannot be described by fluid dynamics. 
We next discuss various theoretical approaches to the entropy growth rate during 
the thermalization stage. 



3. Entropy growth rate 

3.1. Husimi formalism and KS entropy 

If the quark-gluon plasma were a weakly coupled system, the growth of the en- 
tropy during the equilibration phase could be calculated from a partonic Boltzmann 
equation. However, the small shear viscosity and other observations from the RHIC 
experiments, such as measurements of jet quenching, indicate that the plasma is 
strongly coupled. In addition, it is thought that the state created by the deco- 
herence of the nuclear gluon wave functions (the glasma \ \) is characterized by 
strong, still partially coherent color fields. Such fields may even be regenerated by 
plasma instabilities in the pre-equilibrium quark-gluon plasma. A formalism that 
professes to describe entropy production during this early stage must, therefore, 
be able to describe the growing complexity of a quantum system as it makes the 
transition from a regime of gauge field dominance to the hydrodynamical regime 
characterized by complete randomness on thermal length scales. 

A general formalism that can describe this transition in terms of the gauge field 
dynamics has been proposed by Kunihiro et al.^^ The idea underlying this approach 
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is that the growing entropy measures the increasing intrinsic complexity of the 
quantum state of the system after appropriate coarse graining. As mentioned in the 
Introduction, the problem is how to impose a minimal amount of coarse graining 
without assuming the answer. A general solution to this problem dates back to 
Husimi.^^ The Husimi distribution is defined as a convolution of the Wigner function 
X, t) of the system with a minimum-uncertainty Gaussian wave packet: 

H^{p,x;t) = I ^i^expl^-^{p-pr-j{x-xr^W{p',x';t). (13) 

Here x and p stand for all "position" and "momentum" variables characterizing 
the system, which may include particle positions as well as field amplitudes. The 
Husimi distribution depends on the squeezing parameter A, which measures the 
ratio of position and momentum uncertainty. Its value can be chosen at liberty, but 
once fixed, will not change with time. The Husimi distribution can be understood as 
a coarse-grained phase space distribution of the quantum system, where the coarse 
graining corresponds to the projection on a coherent state. We recall that coherent 
states are the closest quantum analogues of classical systems compatible with the 
uncertainty relation. 

Because the Husimi distribution can be shown to be positive (semi-) definite, it 
permits the definition of a coarse-grained entropy, first introduced by Wehrl^ 

SHAt) = -l ^ H a{p, x;t) In H^{p,x;t). (14) 

Quantum systems containing unstable modes, i.e. modes with an exponentially 
growing amplitude, have a linearly growing Husimi- Wehrl entropy.!^ The entropy 
growth rate is found to be independent of the squeezing parameter A and thus in- 
dependent of details of the coarse graining. It is given by the sum of the exponential 
growth rates of all unstable modes. In classical dynamical systems, this quantity 
is known as the Kolmogorov- Sinai entropy, or KS entropy, and defined as the sum 
over all positive Lyapunov exponents of the system: 

k 

The KS entropy is understood to be a measure of the growth rate of the coarse 
grained entropy of a dynamical system starting from a configuration far away from 
equilibrium, after an initial start-up phase during which unstable fiuctuations grow 
to dominance and before it gets too close to its micro-canonical equilibrium.'^ 

3.2. Eigenstate Thermalization 

As stated in the Introduction, the eigenstate thermalization hypothesiJUE (ETH) 
states that an isolated quantum system thermalizes under its own internal dynamics 
if every energy eigenstate contains a thermal component. More precisely, the ETH 
posits^ that the expectation value {"^ eI^I"^ e) of a few-body observable A in an 
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Fig. 5. The Wigner function for the inverted oscillator at t = and t = 2/ A. The horizontal axis 
denotes the scaled position uux; the vertical axis represents the momentum p. 




Fig. 6. Husimi function ( |13| > for the unstable oscillator at t = and t = 2/ A. Note that the extent 
of the distribution in the off-diagonal direction (p — ujx) does not shrink beyond the resolution 
limit set by the Gaussian smearing introduced by the Husimi transform. 

energy eigenstate ^ e of the Hamiltonian of a large, interacting many-body system 
equals the microcanonical average of A at the mean energy E: 

{^e\A\^e) = {A) (16) 

For a system with very many degrees of freedom, the microcanonical average is gen- 
erally agrees with the thermal average for an appropriate choice of the temperature. 
A corrollary of this hypothesis is that a state that is an eigenstate of a few-body 
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observable A will be a superposition of energy eigenstates 

^A = Y.Ce^E (17) 

E 

with probabilities \Ce\'^ that thermally distributed: \Ce\'^ ~ exp{—E/T). As the 
state evolves with time, 

^A{t) = Y,CEe-''^'/^^E (18) 

E 

the phases of the different energy components increasingly diverge from each other. 
The contributions from different energy eigenstates thus effectively decohere at late 
times. Any measurement of a few-body observable on the pure quantum state yields 
a result that is indistinguishable from the results of a measurement on a mixed state 
defined by the diagonal density matri^P^ 

PE,E' = \CEfSE,E' . (19) 

Rigol et a/.f^ argued on the basis of this argument that the coarse grained entropy of 
a pure quantum state immediately after a quench can be calculated as the entropy 
associated with the diagonal energy density matrix in the Hamiltonian after the 
quench: 

S = -Y,\Ce\'H\Ce\') . (20) 

E 

However, this argument does not take into account the time evolution of the phases 
of the different energy components, which only become quasi-random numbers at 
late times. At a time t after the quench, the phases of energy eigenstates within 
an energy band AE{t) ^ h/t differ by less than 7r/2, and thus these components 
of the wave function remain approximately coherent. The time evolution of the 
coarse grained entropy should therefore be approximately given by the following 
expression: 

S{t) = -^\CESH\CEj') (21) 

Ea 

where denotes an energy band of width AE{t) and 



CEAt)= 2^ Ce . (22) 

\E-Eo,\<AE{t)/2 

As time progresses, AE(t) continues to shrink, until it is eventually smaller than 



any spacing between energy eigenvalues of the system. Thus S{t) defined in (21) 



approaches the asymptotic value (20) at late times. However, immediately after the 
quench, S{t) = 0, and the coarse grained entropy only increases gradually as more 
and more energy eigenstate components dephase from each other. We will show in 
the next subsection quantitatively how this mechanism operates in the case of an 
exactly solvable quantum quench. 
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3.3. A toy example: The inverted harmonic oscillator 
3.3.1. Husimi function approach 

The ideas behind the Husimi- Wehrl entropy are best ihustrated by a simple exam- 
ple. We choose the inverted harmonic oscillator because this case is so simple that 
we can analyze it also in a completely different approach built on the eigenstate 
thermalization hypothesis (ETH). For the inverted harmonic oscillator 

n = lp'- \>?x' (23) 

we start from a Gaussian wave packet 

(x|^o>=(|^)''V-^^/^'^ (24) 

and calculate H/\{p,x;t) analytically. The result is shown in Figs, [s] and |6j which 
illustrate how the phase space volume stays constant for the Wigner function but 
increases for the Husimi function because the collapse in the shrinking directions is 
halted by the resolution A. One finds thalP^ 

lim^ = A . (25) 

It is noteworthy that this results is independent of A and thus independent of the 
details of the coarse graining. 



3.3.2. Eigenstate thermalization approach 

Next we show that on gets the same result in the approach based on the eigenstate 
thermalization hypothesis. The basic concept underlying the ETH is that every 
eigenstate of a Hamiltonian contains a thermal state. The initial quantum state is 
a coherent superposition of such eigenstates with the property that this thermal 
component does not contribute to physical observables. If coherence is lost this 
state becomes visible and can dominate at late times ^. For the inverted harmonic 
oscillator we can use the WKB approximation 

1 r^^^/ 



/ 2 cos i dx' pe(x') , , / ^-r , , 

^£; = a/^ ^-^V^ ^ Pe{x) = V2E ^ X^x^ (26) 



For the Gaussian initial state ( 24 ) one expands the initial state in energy eigenstates 
^{x,t) = j dE Ce e-i^* ^e{x) . (27) 
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The time evolution of a narrow energy band /S.E within the initial state is given by 

JE-AE/2 



(27rfiwi;)i/4(27rM|a;|)i/2 

2ri y A xq u) j 

ifi_ilnM_l) 

2h \ A Xq UJ I 



(28) 



At time the phases remain coherent within an energy band of width /S.E = 
h/t. The ETH approach rests on the fact that this band shrinks with time t. The 
coherence is controlled by the right-most sine function which for large times imposes 
the constraint x ~ xqc"^^ and by the term ex.p{iXx'^ /2h) in the exponent. The phase 
coherent bands in the variable x'^ are determined by the exponential, yielding the 
position uncertainty 

A.~A_^ (29) 

For the occupation probability of each coherent energy band this condition implies 
for large enough times t: 

/oo 2 ^-2E/huj 

dx^EAE{x,t) = 2AE , (30) 

-oo V27rhujE 

which implies the normalization dE pE = AE. Note the appearance of the 
thermal weight factor ^-^^jhuj ^ predicted by the ETH. This result can be easily 
extended to include the block structure in x\ 

,rF ._ e(Ax-|x-xoe^*|) e-^^l^ 

with the normalization 

dE dx p{E, x) = l . (32) 
Jo Jo 

From the block matrix of all coherent energy-position bands p{E^ x) it is then 
straightforward to calculate the coarse grained entropy for large times: 

POO POO 

S= dE dx [p{E,x)\np{E,x)] 
Jo Jo 

POO pOO 

^ dE dx [p{E,x)\n{e^^)] = Xt , (33) 
Jo Jo 

thus leading to a constant entropy growth rate dS/dt = A in agreement with the 



results (25) obtained in the Husimi formalism. 

Apparently, dephasing of {AE^ Ax) blocks has a similar effect as phase space 
smearing, at least in this toy example. This suggests that the phase of linear entropy 
growth may be a generic feature of entropy production in quantum field theories. 



12, 2011 0:11 WSPC/INSTRUCTION FILE BigPicturev6 



18 

3.4. Entropy growth rate of non-Abelian gauge fields 

The analysis of the previous sections suggests that entropy production is, in general, 
dominated by a phase of linear growth with time, where the slope is given by the 
Kolmogorov- Sinai' entropy growth rate (in short, the KS entropy). We assume that 
this is also true for highly excited quark-gluon systems like those created in high- 
energy heavy-ion collisions. In this case, the KS entropy can be calculated from the 
dynamics of classical gauge fields. Kunihiro et a/.f^ studied the real-time dynamics of 
classical Yang-Mills fields numerically after replacing continuous space by a three- 
dimensional cubic lattice. On this lattice they investigated the Hamiltonian 

H=^-J2Etixf + \ Yl ' (34) 

x,a,i x,a,i,j 

(x) = d^A^ix) - d.Atix) + J2 r'"A',ix)A-jix) , (35) 

6,c 

where di denotes the central difference operator in the i-direction. The classical 
equations of motion are then 

A1{x) = Et{x) , (36) 
Et{x) = J2 djF^^ix) + E r''Ajix)n(^) ■ (37) 

j b,c,j 

The spectrum of Lyapunov exponents and the Kolmogorov- Sinai' entropy were de- 
termined with two different distance measures. As expected, the slope was found 
to be independent of the distance measure used, as shown in Figj?] 
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In order to convert the results of the lattice simulations to dimensional, physical 
values, one needs to take into account that classical gauge theories have unphysical 
ultraviolet divergencies. This limitation of classical field theory is already encoun- 
tered for point charges in classical electrodynamics. The physical reason is that 
classical point charges do not exist; no particle can be localized better than within 
its de Broglie wave length. Thus the lattice spacing a is not an unphysical parame- 
ter which should be sent to zero as in lattice simulations of the euclidean quantum 
field theory. Rather it determines the length scale below which neglected quantum 
effects become relevant. There are various ways to determine a. For SU(A/'c) gauge 
theories the inverse damping rate for infrared modes is e.g.'^ 

1 ^ 247r ^ 1 

for Nc = 3. Alternatively one can compare the energy density on the lattice and in 
the continuum 

eLattice(T) = 2iN^ - 1) ^ 
6Cont(T)=2(7V2-l) ^T^ 

30 \i/3 1.4 



/ 1.4 



The fact that a should scale like e~^/^ implies that the entropy slope in lattice units 
should scale like e^/^ which it does, as shown in Figjsj 

Kunihiro et al. also confirmed that the results for the Lyapunov exponents 
are volume independent, see Figj9] Finally, in view of the discussion below it is 
worth mentioning that within classical Yang-Mills theory thermalization follows 
the bottom-up rather than the top-down scenario, i. e. soft modes grow stronger 
and thermalize faster than hard modes. More details can be found in Kunihiro et 
al^ The conclusion from these investigations is that the thermalization time in a 
relativistic heavy-ion collision must be approximately 2 fm/c and can hardly be 
shorter than 1 fm/c. 

The success of the classical lattice simulation for the KS entropy of the non- 
abelian gauge field poses the question how generic this behaviour is. The lattice 
regularized Yang-Mills equations describing the dynamics of classical color fields are 
known to be strongly chaotic, and the KS-entropy of the classical Yang-Mills 
field was shown to be a thermodynamically extensive quantity.ED It could be that 
systems with less pronounced ergodic properties show a less clear-cut behaviour. 
Even for quantum chromodynamics there are still many open issues. For example, 
one can extend the lattice approach to include Gaussian fiuctuations.^ It would 
also be interesting to evaluate the increase of the coarse grained entropy associated 
with the growth of the initial quantum fiuctuations around classical gauge fields in 
the colliding nuclei.^ 
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Fig. 9. Time evolution in SU(2) simulation on 4"^, 10^, and 20"^ lattices with the same energy 
density. 
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Finally, let us comment on related work by Berges et a/P^ These authors an- 
alyzed the isotropization of the energy momentum tensor of classical Yang-Mills 
fields. In their Fig. 8, low momentum modes show a behaviour reminiscent of that 
observed by Kunihoro et al. with a similar isotropization time of 1 — 2 fm/c. On 
the other hand, the unexpected behaviour of high momentum modes may not be 
physically significant, because these modes are subject to lattice artefacts. 

4. Entanglement Entropy 
4.1. The basic concept 

The entanglement entropy is defined as follows: Consider a quantum system X 
composed of two complementary subsystems B and B' such that X = B U B' ^ 
which is in a pure state \X). If the state of the subsystem B is not a pure state, 
one says that the quantum states of B and B^ are entangled. The density matrix 
of the subsystem B can be written as pb = Ttb' {\^){^\)^ and the entanglement 
entropy is defined as 

S{B) = -TTB[pB\npB]. (40) 

This definition measures the amount of entanglement of the quantum state of B 
with the quantum state of B' . Obviously, for this definition to make sense S{B) 
must be equal to S{B^) obtained by tracing the density matrix over B. The rela- 
tionship S{B) = S{B^) implies that the entropy can only depend on the area of the 
common surface separating B from B\ which we call A{B) = A{B^). If the entangle- 
ment entropy is proportional to the surface area, as explicit calculations in tractable 
quantum field theories show, this implies that S{B) = hzA{B) = hzA{B^) = S{B^), 
where Hi is the proportionality constant. Entanglement entropy thus contains by 
construction the holographic principle, making it one of the primary objects of 
interest for the thermodynamic properties of event horizons. The entanglement en- 
tropy was first quantitatively studied for a three-dimensional quantum field theory 
by Srednicki, who found that the value of the entropy depends quadratically on the 
ultraviolet momentum cut-off of the field theory.l^ The concept trivially generalizes 
to spatial dimensions other than d = 3. 

The one-dimensional case has been extensively studied for conformal field 
theories. In that case, the entanglement entropy of an interval of length £ in 
the vacuum state of the field theory is given, up to a possible constant, by the 
general formulalM^ 

S{e) = pn{i/a) , (41) 

where c denotes the central charge of the conformal field theory and a is the short- 
distance cut-off. For this simple case, the analogue to the apparent confiict in three 
dimensions between the volume dependence of thermal entropy and surface depen- 
dence of entanglement entropy can be resolved as follows. When the field is not in 
its vacuum state but in thermal equilibrium at T > 0, then S{B) measures not only 
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the entanglement entropy with the surroundings of domain but also the thermal 
entropy contained in B. As expected, the thermal contribution is proportional to 
the volume of B (in one dimension simply the length £ of the interval). The thermal 



generalization of (41) i&lMi£2i 



Sii) = ^ln( — sinhi7reT)], (42) 



which goes over into (41) in the limit <C 1, and approaches the thermal entropy 
^th = {c/3)7rT£ in the limit Ti > 1. 

When the original vacuum state is perturbed by a quantum quench, the entan- 
glement entropy increases with time and eventually reaches the thermal value. In 
the limit of a near critical one-dimensional conformal field theory on an interval of 
length i one 

where tq is the so-called extrapolation length. The final entropy corresponds to the 
temperature T = (4ro)~^ Remarkably, the entanglement entropy reaches its equi- 
librium value after a finite time ^/2, which implies that the entangled information 
travels away from the center of the interval at the speed of light. 

The concept of entanglement entropy has an important application to black 
holes. In the theoretical investigation of the properties of black holes it was real- 
ized early on that thermodynamic consistency requires black holes to have some 
form of entropy, which was suggested by BekensteiiP^l to be proportional to its 
surface 6'bh = ^/4. The relation dE = dM = Tc^S'bh, which holds when some 
object falls into the black hole, requires black holes to have a nonzero temperature, 
e. g. T = I/SttM for a Schwarzschild black hole, which is exactly the Hawking 
temperature. That the latter is of quantum origin (notwithstanding the fact that 
some of its properties can be derived from thermodynamics) underlines the fact 
that entropy cannot be formulated consistently in a classical theory, as usually dis- 
cussed in connection with the Third Law of Thermodynamics. Equating statistical 
entropy and thermodynamic entropy implies that the number of quantum states 
must be proportional to A/ 4. This relation implies that all information contained 
in a black hole is encoded in the area of its horizon, which is a special case of the 
holographic principle ESESI The latter states that the number of degrees of freedom 

understood as the number of independent quantum states describing a region B 
of space-time is bounded by the area A{B) of its boundary.^^ Many detailed inves- 
tigations suggest that this definition may, indeed, avoid all apparent paradoxes of 
information fiow across an event horizon that have been constructed over the years. 

For our discussion, two aspects are worth highlighting: 



• The surface entropy of a horizon can be identified, up to a constant factor, with 
the entanglement entropy. 
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• There exists a fundamental difference between a static black hole and a dynam- 
ically created and evaporating one, as the entanglement entropy (which is often 
identified with the von Neumann entropy) can remain zero for the latter. Thus 
the unitarity of the boundary gauge field theory might find its counterpart in 
in the unitarity of the black hole formation and evaporation process, see e. g. 
Takayanagi and Ugajin.=SSl 

There is a growing consensus that the formation and evaporation of a black 
hole is a unitary process, at least in asymptotically fiat or AdS space-times!^ This 
implies that the von Neumann entropy stays constant throughout the process and 
information about the initial state can, in principle, be recovered from careful mea- 
surements of correlations within the Hawking radiation.^ The apparent (Beken- 
stein) entropy assigned to a black hole is an expression of the inability of a classical 
observer to distinguish between the many micro-states of the true quantum geom- 
etry. In the "fuzzball" pictur^^ of the black hole, the horizon encloses the region 
in which the average micro-state deviates strongly from the classical geometry. The 
finite size of the fuzzball is set by the uncertainty relation, which dictates that a 
micro-state carrying enough energy to represent the mass of the black hole will, in 
general, be spread out in space. One can easily imagine that a phase space smear- 
ing of a micro-state, similar to the Husimi transform of the Wigner function, would 
map the the pure state of the quantum geometry into the entropy carrying, quasi- 
thermal geometry of a classical black hole. The Bekenstein entropy would then have 
to be understood as the coarse grained entropy for a transitory black hole created 
by gravitational collapse as well as the entanglement entropy of an eternal black 
hole. 

Takayanagi and UgajiiP^ have argued that under certain circumstances the 
entanglement entropy of a transient black hole can remain zero. This could happen, 
e.g., if the forming and decaying black hole is microscopic so that it decays into 
a small number of excitations whose detailed phase relationship can be observed. 
They construct a toy model in which "black hole" formation is a periodic process, 
which obviously does not entail entropy growth. This is a special case of the general 
observation that the concept of entropy for a finite system only makes sense on time 
scales much shorter than the Poincare recurrence time. While this is practically true 
for most systems with many degrees of freedom, it does not apply to microscopic 
systems, and it may not apply to simplified systems that have been constructed to 
be mathematically tractable. 

We also note that the von Neumann entropy of an eternal black hole may be 
non-zero. The fundamental difference between a stationary and a transitory black 
hole with respect to their entanglement entropy is illustrated in FigfTo) Before the 
black hole is formed, B and B' cover the full quantum system X = B U B' . Thus 
entanglement entropy can be defined as the minimal surface between B and B' . If 
a black hole develops inside of B' this minimal surface does not change at all (case 
b). As the process is unitary, B' undergoes some non-trivial internal reordering of 
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states, but this does not affect the surface between B and B' . However, in the 
presence of a stationary black hole, X B U B' at early times and one has to 
define entanglement entropy from the minimal surfaces between the three regions 
B, B' and B" (case c). 



B 



B' 



B 



B' 



B 



B' 



a) 



b) 



Fig. 10. An illustration of the meaning of unitarity within the formulation of entanglement entropy: 
a) two regions of space separated by a boundary; b) one region (B^) containing a transitory black 
hole; c) one region (B^) containing an eternal black hole, which can be considered as a third region 
of space (B^^) with its event horizon as natural boundary with respect to region B^ 



4.2. The AdS/CFT Duality 

The reason we discussed the question of entropy of a black hole at some length is 
that black holes play an important role in holographic gravity duals of interacting 
quantum field theories. The AdS/CFT dualit'^^i^ posits that superstring theory in 
AdSs X is holographically dual to a conformal field theory in Minkowski space, the 
Af = 4: supersymmetric SU {Nc) gauge theory. In the limit of a large number of colors 
Nc the string theory approaches the classical limit of the supergravity theory on 
AdS space. When a black hole with Schwarzschild radius much larger than the AdS 
curvature radius (a "black brane") is added to the AdSs space, the string theory 
becomes dual to the thermal super- Yang- Mills theory in Minkowski space. The 
connection to relativistic heavy-ion collisions originates in the observation that for 
very high temperatures, i. e. deeply in the deconfined phase, QCD is approximately 
conformal and fermionic contributions to thermal properties should be much less 
important than gluonic ones for large Nc. In addition, supersymmetry is broken by 
finite temperature, and the adjoint superpartners of the gauge field ("gauginos") 
become effectively massive and decouple from the infrared sector of the theory. Thus 
the absence of A/'-fold supersymmetry in thermal QCD should be less relevant than 
in the vacuum theory. Finally, the question whether the thermodynamic properties 
of SU{Nc) gauge theories already have reached the large- A/'c limit for Nc = 3 has 
recently been answered in the affirmative by lattice-QCD simulations.^^ 

For those readers who are not familiar with AdS/CFT duality, we now give a 
very brief description and refer to Natsuume^^ for more details and references. As 
already stated, AdS/CFT duality establishes a correspondence between superstring 
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theory in the AdSs x space-time and the A/" = 4 supersymmetric SU (Nc) gauge 
theory in Minkowski space. For many purposes, the part of the ten-dimensional 
space-time can be ignored. The correspondence is thus between a five-dimensional 
theory of quantum gravity and a four-dimensional gauge theory. The most im- 
portant aspect of the correspondence is that the strong 't Hooft coupling limit 
{g'^Nc oo) of the gauge theory is mapped onto the weak coupling limit of the 
string theory, which is well described by its classical limit, the supergravity theory 
on space-times with asymptotic AdSs geometry. The field theory can be consid- 
ered as "living" on the asymptotic boundary of the AdSs space-time; the AdSs 
space-time is therefore usually referred to as the "bulk" . 

There is an extensive dictionary of correspondences between expectation values 
of operators in the field theory and geometrically defined objects in the dual super- 
gravity theory.l^For example, the expectation values of some local operators of the 
gauge theory, such as the energy-momentum tensor, are related to the metric or 
other local fields in the gravity theory. The expectation values of non-local opera- 
tors in the gauge theory, such as two-point functions or Wilson loops, correspond to 
the actions associated with geodesies, minimal surfaces, or other geometric objects 
in the dual supergravity theory. 

The thermal state of the field theory, i. e. the equilibrated gauge field plasma, 
is represented by Anti-deSitter space with an imbedded black hole of mass M, 
whose Hawking temperature Th — M^l^ jn equals the physical temperature T of 
the field theory. Empty Anti-deSitter space corresponds to the vacuum state of the 
field theory. Thermalization of the theory is thus dual to the process of black hole 
formation in the bulk. The energy deposition occurring in a heavy ion collision 
corresponds to the injection of energy in the outer regions of the bulk, which then 
falls by gravitational attraction and eventually forms a black hole, providing a 
holographic description of the process of thermalization in the gauge theory. 

In recent years it has been understood that AdS/CFT duality makes it possible 
to obtain rigorous predictions for strongly coupled gauge theories at high tempera- 
ture which seem to agree with phenomenological results deduced from high-energy 
heavy- ion collision data.^^ However, it is important to keep in mind that QCD is 
not a strongly coupled theory in the high temperature limit, where it becomes ap- 
proximately conformal. On the other hand, in the temperature range of practical 
interest (T < 3Tc ~ 500 MeV) where thermal QCD is a strongly coupled theory, 
it is far from being conformally invariant. Lattice-QCD simulations show that the 
trace anomaly of the energy-momentum tensor, = e — 3P, is large and of the 
same order as e itself (see Borsanyi et al.^^^ Fig. 6). It is thus not obvious that the 
AdS/CFT correspondence in its simplest form is a good model for thermal QCD in 
the region of interest. This has motivated attempts to find improved holographic 
duals (see, e. g., Alanen et a/.'^, Kajantie et al.^^). 
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4.3. Holographic Entanglement Entropy 

A volume of Minkowski space V filled with a quantum field in its vacuum state has 
a non- vanishing entanglement entropy associated with it, because the field modes 



of its complement necessarily leak across the boundary (see Fig. 11). Although the 
field is known to be in the vacuum state within V, nothing is assumed to be known 
about the state of the field outside. As mentioned in Section |4.1[ this ultraviolet 
divergent entanglement entropy is proportional to the surface area of the volume.^^ 
A holographic dual description of the entanglement entropy of a region of Minkowski 
space-time was proposed by Ryu and Takayanagi^^ (see Nishioka et al.^^ for a 
review). They extended the surface of the volume into dimension of the bulk to 
form a minimal hypersurface 7(V^), whose boundary coincides with the surface of 
the Minkowski space volume V. This is illustrated in Fig. [TT] The entanglement 
entropy associated with this surface is defined asHHHl 

- -4G^ ^^^^ 

where ||7(^)|| denotes the volume of 7(V), which has the same dimensionality as 
V, and Gn is Newton's gravitational constant. When the field is in the vacuum 
state, j{V) extends into pure AdSs space, and one finds that its ultraviolet diver- 
gent volume is proportional to the surface area of V. When the Anti-deSitter space 
contains a black hole and the region V is large enough, the hypersurface j{V) hugs 
the event horizon of the black hole, and ||7(V^)|| acquires an additional contribution 
proportional to the covered horizon area, which is proportinal to the volume of V. 
By explicit calculation one finds that the additional contribution to the entangle- 
ment entropy is precisely the thermal entropy of the gauge field in the Minkowski 
space region V. Ryu an Takayanagi confirmed the equivalence between the entan- 
glement entropy calculated in the field theory and that calculated holographically 



by means of (44) in several other tractable cases, but no explicit equivalence proof 



is known. In the following, we assume that the equivalence holds, in general. 



The relation (44) is remarkable in that it relates the Bekenstein entropy of black 
holes to the usual thermal entropy of a field theory by holographic duality. We can 
imagine applying this concept to a relativistic heavy ion collision, in which a certain 
region of space is suddenly filled with high energy that quickly equilibrates into a hot 
quark-gluon plasma, which then expands, cools, and eventually dissolves. In the dual 
description, energy is injected asymptotically into an AdSs space, collapses to form 
a black hole (or black brane), which eventually evaporates by Hawking radiation. 
The formation of the black hole, i. e. the dual of the thermalization process, can 



be tracked by the growth of the entanglement entropy defined in (44). However, in 
this context, the holographic duality of the entanglement entropy concept raises a 
number of questions: 



• How is the entropy of a black brane in AdS^ related to the von Neumann and 
Wehrl entropies in the boundary gauge theory? 
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surface 



Fig. 11. Holographic representation of the entanglement entropy of a finite region V of space as 
the area of the minimal hypersurface j{V) in the bulk, which has the same boundary as V. 

• Can the duality be extended to the dynamical process of black brane formation 
in AdSs and thermalization in the boundary field theory? 

• Is it correct to apply the entanglement entropy correspondence, which has been 
verified to be dual in the case of thermal equilibrium, to non-equilibrium situ- 
ations? 

None of these questions can be conclusively answered at present. Nevertheless, we 
want to review some of the ideas and insights which might contribute to find relevant 
answers in the future. Also we will review results from one specific model study. 

A crucial element in the discussion of thermalization is the form the fluctuation- 
dissipation theorem takes in the bulk. Following Caron-Huot et alJ^ one can argue 
as follows: The large black brane in AdS emits Hawking radiation in the same way 
as black holes in Minkowski space. This does not lead to a heating of the boundary, 
because perturbations along outgoing geodesies are reflected at the boundary and, 
therefore, end up again in the black brane. However, the quantum fluctuations 
associated with this Hawking radiation induce fluctuations in the boundary, which 
by the fluctation-dissipation theorem are related to dissipation, and push the system 
towards equilibrium. The closer a fluctuation is located to the horizon, the longer 
it takes to reach the boundary and the stronger it will be red-shifted, while all 
higher momentum fluctuations get reabsorbed by the black brane on a time scale 
of order 1/T. In the boundary field theory, therefore, ultraviolet modes thermalize 
first, infrared modes last. 
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To gain more insight into this complex situation one studies the reaction of 
AdS/CFT for various types of sudden violent perturbations, often called quan- 
tum quenches. The entanglement entropy then traces the evolution from such a 
highly excited, but highly phase correlated initial state to the state of thermal 
equilibrium!^ To reduce the technical difficulty of such calculations, one often stud- 
ies lower dimensional theories. In the next section we will review results in 2, 3, 
and 4 dimensions for a schematic model of black hole formation in the bi 

4.4. Thermalization 

Albash et al7'^ and Balasubramanian et a/P^Hl considered the case where a finite 
energy density was deposited in the bulk near the boundary at a certain initial 
time, resulting in a shell of null dust falling in the fifth dimension until it forms an 
event horizon and a black brane is created. The resulting metric was then probed 
by strings and membranes with endpoints on the boundary. Thermalization was 
studied as a function of the Minkowski space separation between the endpoints, 
respectively the area of the Wilson loop, as a function of time. The time required 
for the geodesic length or the membrane area to reach that of the black brane equi- 
librium configuration was interpreted as an estimate for the physical thermalization 
time, as measured by the entanglement entropy. The results of this study showed 
that complete thermalization occurs within a time that is equal to half of the spatial 
diameter of the probe in the boundary space. 

The AdS space dual to the thermalizing field theory is probed by geodesies 
and surfaces with endpoints on the boundary, see FigjT2] The effect of the in-falling 
shell can be expressed geometrically in so-called Poincare coordinates by the Vaidya 
metric 

ds^ = ^ [- (1 - rn{v)z'^) dv'^ - 2dzdv + dx^] , (45) 

where v labels in-going null trajectories and the AdS radius has been set equal 
to unity. The boundary space-time is located at z = and x = {xi, . . . , Xd-i) 
correspond to the spatial coordinates on the boundary. The mass function of the 
in-falling shell was taken to be: 

M / v\ 
m{v) = ^ 1 + tanh — , (46) 

where vq parametrizes the thickness of the shell falling along v = 0, which is assumed 
to be very small E3Z§EllZl] For a sufficiently large separation of the two endpoints on 
the boundary and before the shell has formed an event horizon, the geodesic curve 
in the bulk connecting the endpoints on the AdS boundary will "punch" through 
the shell and explore the still un-thermalized AdS geometry below the shell. The 
larger the portion of the curve extending below the shell, the larger is the difference 
between the geodesic length of the curve and that of a geodesic connecting the two 
endpoints in the AdS-Schwarzschild geometry, which represents full thermalization. 
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Fig. 12. Illustration of the case studied. Energy is deposited on the boundary at a certain starting 
time. This situation should be seen as model for the creation of a fireball in a high-energy heavy 
ion collision. This energy deposition leads to the propagation of a shell of what is called 'null 
dust' in the fifth dimension. Its dynamics can be solved analytically and absorbed into a specific 
metric, the Vaidya metric. The latter is then probed with e.g. a string with fixed endpoints on 
the boundary, which corresponds to the correlation function between two operators of very large 
mass-dimension. The analytic form of the geodetic is known for the equilibrium case after the 
shell moved beyond the horizon. The time scale with which the geodetic length approaches its 
value for that solution is identified with the thermalization time. 




Fig. 13. Example of a space-like geodesic that starts and ends on the boundary of AdS (z = 0) 
with a separation xq. Outside the shell, the geodesic propagates in a black brane geometry, while 
inside it propagates in an empty AdS geometry. 



It is then possible to obtain analytic or numerical solutions for different dimen- 
sions and probes. Some results were obtained by Balasubramanian et a/ElEll for 
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the entanglement entropy. In one spatial dimension {d = 2) the entanglement en- 
tropy is holographically given by the length of the bulk geodesic connecting the 
two endpoints of the linear region in the boundary space. In two spatial dimen- 
sions {d = 3) the entanglement entropy is holographically given by the minimal 
area of the bulk surface, whose perimeter coincides with the perimeter of the cir- 
cular boundary area. In three spatial dimensions {d = 4) the entanglement entropy 
is holographically given by the minimal volume of the bulk hyper surf ace, whose 
surface coincides with the surface of the spherical boundary region. 

The results are shown in Fig. [l4j In each dimension, the difference between 
the minimal length (or area or volume) in the bulk for the geometry at a given 
time and minimal length (area, volume) in the fully thermalized AdS-Schwarzschild 
geometry is shown, normalized to the boundary volume {jC = C/£^ A = A/((nR^)^ 
V = V/(47ri?^/3)). In every case one finds that thermalization as measured by 
the entanglement entropy is achieved within the time given by precisely half the 
diameter of the probe in the boundary divided by the speed of light. The main 
results of this investigation arelI3I3 

• Thermalization is an extremely fast process in the strongly coupled gauge the- 
ory, which is only constrained by causality. The information loss from a finite 
region proceeds at the speed of light. This is especially remarkable because the 
strongly coupled gauge theory does not support quasi-particle excitations, with 
the exception of phonons, which do not play a role in this process. The mecha- 
nism by which the information is transported out of the thermalizing volume is 
unknown. However, the result suggests that the strongly coupled gauge theory 
behaves like a fast scrambler as defined by Sekino and Susskind.^^ 

• Other observables reach their equilibrium values faster than the entanglement 
entropy. This is not surprising, because the entropy is sensitive to all degrees 
of freedom of the gauge theory, while other observables are sensitive only to a 
subset. 

• Short distances thermalize first, i.e. thermalization occurs top-down for a gauge 
theory in the limit of very strong coupling, while it occurs bottom-up for weakly 
coupled gauge theories.^ 

5. Summary 

The problem of thermalization of highly excited states in quantum chromodynamics 
and gravity has been intensively studied during the past decade. In QCD, the 
interest was driven by the experimental results from relativist ic heavy ion collisions, 
which indicated that a thermal quark-gluon plasma is formed on an extremely 
short time scale of order 1 fm/c. In gravity theories, the main motivation was the 
desire to resolve the information paradox associated with black hole formation. In 
both situations, the question can be phrased in terms of the growth rate of an 
appropriately defined entropy. While the von Neumann entropy of a closed system 
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Fig. 14. The difference in normalized geodesic length A(6jC), minimal surface area A(6A), and 
minimal volume A(6V), respectively, to that of the thermalized state as a function of boundary 
time to for dimensions d = 2 (left), d = 3 (center), and d = A (right). The results are for a thin 
in-falling shell {vq = 0.01). The boundary separations were taken to be ^ = 1,2,3,4 and radii 
R = 0.5, 1, 1.5, 2 (from top to bottom curve), respectively. All quantities are given in units of M. 
Complete thermalization occurs precisely at = ^/2 or Iq = R, respectively. 



does not change with time in any quantum theory due to the unitarity of the 
time evolution, other measures of entropy, which take into account the principal or 
practical inability of an observer to measure all details of a system, can grow with 
time and eventually approach the thermal equilibrium value. 

In this review, we covered several different definitions of such a coarse grained 
entropy with broad applicability: Husimi's phase-space smearing of the Wigner dis- 
tribution, de-phasing of the quantum state underlying the eigenstate thermalization 
hypothesis, and entanglement entropy applied to holographic duals of the quantum 
system. We showed that de-phasing and Husimi smearing result in the same en- 
tropy growth rate in the case of a single unstable mode, which was equal to the KS 
entropy describing the growth rate of the coarse grained entropy in the classical 
limit. This observation motivated a numerical study of the entropy growth rate of 
the non-abelian SU(2) gauge theory, reviewed in Section 3.4, which confirmed the 
rapid thermalization observed in the experiments. 

The concept of entanglement entropy and its holographic equivalent enables 
rigorous calculations of the approach to thermalization in the strongly coupled 
conformal supersymmetric gauge theory. In this approach, the process of thermal- 
ization is dual to the formation of a black hole, linking the study of thermaliza- 
tion in the (3 + l)-dimensional gauge theory to the study of black hole formation 
in (4 + l)-dimensional supergravity. The concept of entanglement entropy for a 
bounded region of space parallels that of the Husimi smeared entropy (the Wehrl 
entropy), because it accounts for the position uncertainty of field modes with a 
given momentum. The two concepts differ in that the entanglement entropy only 
accounts for the uncertainty relation smearing of field modes near the boundary, 
while the Wehrl entropy implements the uncertainty smearing of all field modes. 
The holographic duality also provides compelling evidence for the view that the 
formation and subsequent evaporation of a black hole must be a unitary process 
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in a quantum theory of gravity, which does not result in an increase of the von 
Neumann entropy, but is accompanied by a maximal increase in an appropriately 
defined coarse grained entropy. 

Notwithstanding these significant theoretical advances, there remain many un- 
resolved questions. How exactly are the different notions of coarse grained entropy 
connected? When does black hole formation and decay result in an increase in the 
coarse grained entropy, and when does it represent a practically reversible process 
as in the toy model constructed by Takayanagi and Ugajin?^^ How do the same 
concepts apply to the cosmic Big Bang? We hope that this review will motivate 
some of its readers to study and ultimately resolve such questions. 
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